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Abstract 

Modeling of ion-implanted boron redistribution in silicon crystals 
during low-temperature annealing with a small thermal budget has 
been carried out. It was shown that formation of "tails" in the low- 
concentration region of impurity profiles occurs due to the long-range 
migration of boron interstitials. 
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1 Introduction 

It is well known that device scaling requires decreasing all vertical and lateral 
dimensions of transistors At present, a low-energy high-dose ion implanta- 
tion is used for the formation of shallow junctions However, we need to 
undergo silicon substrate to subsequent annealing to avoid defects created 
by ion implantation. During thermal treatment of ion implanted layers a 
great amount of nonequilibrium self-interstitials is generated, which results 
in transient enhanced diffusion (TED) of impurity atoms. Intensive TED pre- 
vents the formation of shallow junctions especially for boron doped regions 
(see, for example, [U El El II])- Therefore, the problem of formation of doped 
regions with abrupt impurity distributions in vertical and lateral directions 
requires annealing imposed with an extremely small thermal budget for the 
TED to be decreased j5]. One of the techniques to provide a small thermal 
budget and solve this problem is the use of low-temperature annealing of 
short duration [6]. Another technique widely used to suppress transient en- 
hanced diffusion of ion-implanted impurity is the introduction of boron into 
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the silicon layer preamorphized by implantation of Ge ions [TJ, [HI El EE] • As a 
result of the solid phase epitaxial recrystallization (SPER) of the amorphous 
layer at the initial stage of annealing a perfect crystalline structure without 
defects visible by means of electronic microscopy is formed in the array of the 
implanted boron. Nevertheless, in this case TED is also observed during the 
subsequent annealing, but of absolutely another character, and characterized 
by a smaller intensity. Indeed, at low annealing temperatures an extended 
"tail" in the low-concentration region of boron profile is formed in all the 
cases mentioned above. This "tail" represents a direct line for logarithmic 
scale on the concentration axis. 

It was shown in [HJ[T2J[T3] that the majority of the cases of "tail" forma- 
tion in the low-concentration region of ion-implanted impurity distribution, 
including boron implantation in preamorphized silicon, is related to the phe- 
nomenon of a long-range migration of nonequilibrium impurity interstitials, 
and we can neglect channeling of a fraction of implanted atoms if ion implan- 
tation of In, P, and Ga is performed in a random direction or B is implanted 
in a preamorphized layer. In [HI [TJ] , a similar theoretical investigation and 
simulation of impurity redistribution were carried out for boron implantation 
in a random direction in the crystalline silicon. It has been shown that the 
long-range migration of nonequilibrium boron interstitials is responsible for 
the "tail" formation in a crystal lattice damaged by boron ions implanted at 
medium fluencies which are insufficient to create an amorphous layer. For 
example, it was found in [T3] that the average migration length of boron 
interstitials Iai is equal to 24 nm for annealing at a temperature of 600 °C 
for 10 s. This value is two times greater than the average migration length 
of boron interstitials Iai = 12 nm and Iai = H rim found for boron diffu- 
sion in the preamorphized and recrystallized silicon layer at temperatures 
of 850 °C [ID and 800 °C d Ej with annealing duration of 60 s. It is 
interesting to note that due to the higher temperatures and longer duration 
a significant increase in the thermal budget occurred in the experiments that 
were simulated in [11] and [T6| ITT] . Therefore, it is reasonable to investi- 
gate the mechanism of "tail" formation for low-energy boron implantation 
in crystalline silicon during the subsequent low-temperature annealing with 
durations longer than 10 s, i. e., for the increased thermal budget. It is 
worth noting that the elucidation of the mechanisms of "tail" formation has 
a practical importance because this phenomenon influences the depth of p-n 
junction, as well as it is actual needed for the investigation of the transport 
processes of implanted boron atoms. Indeed, it has been shown in [T6"| [T7] 
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that the mechanism of boron diffusion varies with a change in the annealing 
temperature. 



2 Model of boron diffusion 

For simulation of interstitial diffusion of boron atoms, we use the model 
proposed in [Hi [12], which includes the following system of equations: 

1) the conservation law for immobile impurity atoms 

dC T (x,t) _ C AI (x,t) AI 

Q t ~ T AI ^ V X '''V ' \ L ) 

2) the diffusion equation for nonequilibrium boron interstitials 

^^--^u+^,0 = , (2) 
which can be present in the form suitable for numerical or analytical solution 

IB], 

d 2 C AI C AI ~ g AI (Xjt) 

"xr - j2~ + —72 - u > y 6 ) 

OX l M L AI 

where 

l AI = Vd^Ai, ~g AI (x,t)=G AI (x,t)r AI . (4) 

Here C T is the total concentration of substitutionally dissolved impurity 
atoms, impurity atoms incorporated into clusters and trapped by various 
defects including the extended defects (immobile boron atoms); C AI is the 
summarized concentration of nonequilibrium boron interstitials in different 
charge states; d AI and t ai are the diffusivity and average lifetime of nonequi- 
librium impurity interstitials, respectively; Iai and G AI are the average mi- 
gration length and generation rate of interstitial impurity atoms, respectively. 
The concentration of impurity interstitials (mobile species) is also included 
in C T . 

We use the stationary diffusion equation for interstitial impurity atoms in 
view of their large average migration length l AI {l AI ^> lf a u) and due to the 
small average lifetime of these nonequilibrium interstitials t ai {t ai ^ r. 



pr ) 



Here lf a u is the characteristic length relevant to the significant decrease of im- 
purity concentration in the high-concentration region and r pr is the duration 
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of thermal treatment. It is also supposed that the function G AI (x,t) that 
describes the generation rate of interstitial impurity atoms changes slowly 
with time, increasing as compared to the value of r AI . 

It is worth noting that according to the experimental data [TJ [19], a 
significant part of boron atoms in the "tail" region is electrically inactive, i.e., 
impurity atoms do not occupy the substitutional position. A few possible 
explanations can be given to this phenomenon. For example, the average 
lifetime of the nonequilibrium boron interstitials r ~ r pr , i.e., it is equal 
to tens of minutes at temperatures of 800 °C p] and 900 °C [19J . However, 
this assumption contradicts the common point of view that the interstitial 
impurity atoms are highly mobile Really, Iai is not too large (Iai ~20 nm), 
whereas the duration of annealing is longer than enough and for the high 
value of diffusivity of boron interstitials a significant value of Iai should be 
expected. In our opinion, the assumption that the boron interstitials in the 
"tail" region are being incorporated with the defects of the crystal lattice 
and, therefore, become immobile is more reasonable. However, this problem 
requires further investigation. 

In contrast to the models of [TTj [T2] . the Pearson type IV distribution 
[20] is used to describe the spatial distribution of impurity atoms after im- 
plantation: 



where 



C (x) = C{x, 0) = C m F p (x, R p , AR P , Sk) , (5) 



C m = jS— xlO- 8 [^m- 3 ]. (6) 



v^nARp 

Here F p (x, R p , AR P , Sk) is the function that describes the Pearson type 
IV distribution and which is tabulated in [20]; Q is the dose of ion implan- 
tation [ion/cm 2 ]; R p and AR p are the average projective range of implanted 
ions and straggling of the projective range, respectively; Sk is the skewness 
of the Pearson type IV distribution 

On the other hand, to simplify the description for the spatial distribution 
of generation rate of nonequilibrium boron interstitials we used the Gaussian 
distribution: 



G AI (x,t) = G£ 7 (t)exp 



\X Rm] 

2AR* 



(7) 
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shifted to the position of R m to adjust with the maximum of distribution of 
implanted ions. Here R m is the position of a maximum of the Pearson type 
IV distribution for impurity concentration; G^(t) is the value of the gen- 
eration rate for interstitial impurity atoms in the maximum of the Pearson 
type IV distribution. As well as in the models of jTH [12] , it is supposed that 
boron atoms occur in the interstitial positions either directly during implan- 
tation, or during rearrangement or annealing of the clusters that incorporate 
impurity atoms. 



3 Analytical solution of diffusion equation 

It can be seen from expression (j7J) that the generation rate of nonequilibrium 
boron interstitials is approximated by the product of two multipliers, one of 
them describes the dependence of the generation rate on time, whereas the 
other describes the spatial distribution of the generation rate of nonequilib- 
rium boron interstitials: 

G AI (x,t)=Gi I (t)Gf<(x), (8) 

Let us substitute expression ([8]) into Eq.(T5]) and integrate the resulting 
equation with respect to time over the interval [to,tk]: 



ri 2 r AI 

d AI d U eff 

dx 2 



U eff 



T 



AI 



<-fc 

(t k -t ) + Gf*(x)j 'c A J(t)dt 



0, (9) 



to 



or 



AinAi pAi 

d AI -^-^f + 9^(x) = 0, (10) 



dx 2 T 

where 



C eff( X ) 



(tk — to) 



C AI (x,t) dt. 



to 

tk 



o AI 

■J rnk 



(tk — 1 ) 



(11) 



(12) 
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Here C^L(x) and g^ k are, respectively, the concentration of impurity 
interstitials and the value of the generation rate for boron interstitials in 
the maximum of space distribution which are time-averaged on the interval 

[tojtk]- 

It follows from the transformations performed that for the case of the 
generation rate for boron interstitials, representing the product of time- and 
space dependent multipliers, the system of Eqs. (jlj and (J2J) is similar to the 
system of equations (pQ) and (TTOl) with the constant coefficient g AI k . It is worth 
noting that Eq. fllO]) can be solved analytically. For example, such solution 
was obtained in [18] for 



AI r*AI, 



sp 



n A1 r 

9mk ^ 



X 



A T 

9 mk ex P 



x 



2AR 2 p 



(13) 



It is well known that the intensity of transient enhanced diffusion de- 
creases with time. Let us suppose that 



G m \t) = GZ exp 



AI 



t-t 



'gen 



(14) 



where r gen is the relaxation time for the generation of nonequilibrium impu- 
rity interstitials and G^ is the value of the generation rate in the maximum 
of the Pearson type IV distribution for t = to- 

The parameters G^o an d r gen can be derived if we obtain two values of g^ k 
for two thermal treatments with different durations. Indeed, let us suppose 
for simplicity that to =0 and and g^ 2 are the values of the generation rate 
for boron interstitials in the maximum of space distribution time-averaged 
on the intervals [0,ti] and [0,t 2 ], respectively. Then, the following system of 
two independent algebraic equations is valid: 



,AI 



yAI 



9ml ^mO T gen 



1 — exp 



T, 



gen 



(15) 



a AI 

9m2 



r A1 t 

"mO 'gen 



exp 



r, 



gen / J 



(16) 



To obtain the solution of the system fTl5|) . f|T6l) . let us specify the quantity 



kasA = g% x t x /(g A1 
the nonlinear equation: 



,a2 ■ 



Then, the value of r gen can be obtained by solving 
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-^-)-Aexp(-^-)+(A-l) = . (17) 

Miiller's iteration method of successive bisections and inverse parabolic 
interpolation [21] was used in this work to solve Eq. fjTT]) . Then, we obtained 
the value of G^ from Eq. (fT5|) or Eq.flTSD. 

4 Simulation results 

The results of simulation of boron interstitial diffusion within the framework 
of the model described above are presented in Figs. [1] and [2j For comparison 
we use the experimental data of [6J. In paper [6], implantation of boron 
ions was carried out with a dose of lxl0 14 cm -2 and an energy of 1 keV in 
silicon substrates of orientation (100) under 7° tilt, 30° twist angles to prevent 
channeling. Before implantation a pure silicon layer (10-/zm thick, p-type 
conductivity with a resistivity of 10 Qm) has been grown up epitaxially. The 
chemical profiles of the samples before and after annealing were measured by 
means of the high depth resolution secondary ion mass spectrometry (SIMS). 
The annealing of the samples was carried out at a temperature of 600 °C for 
10, 600, and 1200 s. 

It follows from the experimental data of [6] that the boron concentra- 
tion profile after implantation is characterized by a significant asymmetry. 
For this reason, the Pearson type IV distribution was used to calculate the 
concentration profile of ion-implanted boron. It was found in [T5] that the 
values of R p = 0.007 jum, AR P = 0.0055 //m, and Sk = 1.2 give the best 
fit to the experimental data of [6]. It is necessary to note that the Pearson 
type IV distribution with obtained values of the parameters gives a good fit 
for the whole impurity concentration profile excluding the low-concentration 
region. The problem of the accurate calculation of boron concentration pro- 
file in the low-concentration region has been also considered in [15] where 
it was supposed that heating of silicon substrate takes place during the ion 
implantation described in [6]. Then, due to the elevated temperature, a fast 
interstitial diffusion of boron atoms that occupy interstitial position occurs 
which results in the formation of the "tail" in the low-concentration region 
of boron profile. It was found in JTH] that the following values of the pa- 
rameters describing the process of interstitial diffusion allow us to obtain a 
good fit to the boron profile in the low-concentration region, namely: the 
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Figure 1: Calculated boron concentration profile after annealing for 600 s at 
a temperature of 600 °C. The dotted curve represents boron concentration 
profile after ion implantation at an energy of 1 keV with a dose of lxlO 14 
cm -2 described by the Pearson type IV distribution. The dashed curve rep- 
resents boron concentration profile calculated under the assumption of the 
long-range interstitial migration of a part of the impurity atoms during ion 
implantation to provide fitting with the experimental data for boron distri- 
bution after implantation [6] . The experimental data for boron concentration 
after annealing (filled circles) are also taken from [6] 

average migration length of boron interstitials Iai = 0.024 fim (24 nm) and 
a part of boron atoms participating in the interstitial diffusion during im- 
plantation pai is approximately equal to 0.3 %. In addition, it was shown 
in [15] that during subsequent annealing at a temperature of 600 °C for 10 
s the interstitial boron diffusion is characterized by the same value of the 
average migration length of impurity interstitials Iai — 24 nm. The best fit 
to the boron concentration profile after annealing was obtained under the 
assumption that approximately 0.7 % of the implanted boron atoms entering 
into the interstitial position in the high-concentration region of ion-implanted 
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Figure 2: Calculated boron concentration profile after annealing for 600 s at 
a temperature of 600 °C. The dotted curve represents boron concentration 
profile after ion implantation at an energy of 1 keV with a dose of lxlO 14 
cm -2 described by the Pearson type IV distribution. The dashed curve rep- 
resents boron concentration profile calculated under the assumption of the 
long-range interstitial migration of a part of the impurity atoms during ion 
implantation. The experimental data for boron concentration after annealing 
(filled circles) are taken from [6] 



impurity participated in the fast interstitial diffusion during annealing and 
then became immobile again, occupying substitutional positions or forming 
complexes with the crystal lattice defects. 

Taking into account the results of [15], in this work we used the values of 
R p = 0.007 /im, AR P = 0.0055 /im, Sk = 1.2, l AI = 0.024 /im, and p AI = 
0.308 % to calculate the boron concentration profile after implantation. For 
this value of skewness Sk, R m = 0.00521 /im. As can be seen from Figs. [Hand 
121 the results of simulation of boron diffusion during annealing for 600 and 
1200 s agree well with the experimental data if the same value of the average 
migration length of impurity interstitials Iai = 24 nm is used. Taking into 
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account the results of [TJ)] and the present calculations, one can conclude that 
the migration length of impurity interstitials varies negligibly with increase of 
annealing duration. A good fit to the experimental data obtained for different 
annealing durations with the same value of Iai confirms the correctness of 
the developed model of the long-range migration for boron interstitial atoms 
in ion-implanted layers during low-temperature annealing. 

5 Discussion 

It is worth noting that in paper [22] the error arising due to the use of Gaus- 
sian distribution instead of the Pearson type IV distribution when modeling 
the redistribution of ion-implanted impurity was investigated. It was sup- 
posed that the Fick's law describes the diffusion of impurity atoms. As 
follows from the results of [22J , the error in calculating the impurity concen- 
tration profile is less than 5% if L > 5\R m — R p \. Here L = a/ D r pr C AI / Cf 1 
is the average length of the boron redistribution due to the pair diffusion 
mechanism; D is the diffusivity of impurity atoms; C AI is the concentration 
of nonequilibrium silicon self-interstitials; Cf 1 is the thermally equilibrium 
value of self-interstitials. Extrapolating the results of [22J to the case of 
diffusion of impurity interstitials and taking into account that for the exper- 
imental data of [6] simulated in this paper \R m — R p \ is equal to 1.79 nm 
whereas Iai = 24 nm, we can conclude that using the Gaussian distribution 
to describe the generation rate of boron interstitials is quite reasonable. In 
addition, in our computations the Gaussian distribution of generation rate 
is shifted to the position of R m and, therefore, the error in calculating the 
concentration profile for nonequilibrium boron interstitials is quite negligible. 

The best fit to the boron concentration profiles after annealing presented 
in Figs. [T]and [2] was obtained under the assumption that approximately 1.5 
% and 2.29 % of the implanted boron atoms participated temporally in the 
fast interstitial diffusion during annealing for 600 and 1200 s, respectively. It 
follows from these values and the value of pai = 0.7 % for the part of boron 
atoms participating in interstitial diffusion during annealing for 10 s [15] that 
there is a decrease in the generation rate of boron interstitials with increase 
in annealing duration. Indeed, the time-average interstitial generation rates 
for the intervals [0,10] s, [0,600] s, and [0,1200] s are g£ 10 = 6.3xl0 5 , g^ 600 
= 2.19x10 s , and 

5^1200 — 1.675xl0 3 jum 3 s 1 , respectively. The solution of 
the system (TIB]) . (JIB)) for the time intervals [0,10] s and [0,600] s gives the 
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following values of the parameters of time-dependent function ( 1141) for the 
generation rate of boron interstitials: G^ = 8.5655xl0 4 /!m _3 s _1 and r gen = 
15.33 s. The kinetics of generation rate for nonequilibrium boron interstitials 
calculated for these values of G^ and r gen is presented in Fig. [3j A similar 
calculation for the intervals [0-10] s and [0-1200] s gives the following values 
of the kinetics parameters: G^o = 7.5473xl0 4 /xm" 3 s _1 and r gen = 26.6 s. It 
follows from the values obtained that the main fraction of boron interstitials 
is generated during annealing duration shorter than 30 s. Therefore, the 
additional experiments with annealing durations in the range of 10 - 600 s 
are needed to improve the accuracy of kinetics determination. 




Figure 3: Calculated kinetics of the generation rate of interstitial boron atoms 
during annealing at a temperature of 600 °C. The filled circle represents the 
time-average generation rate of boron interstitials for the interval from 
s to 10 s. The dashed curve with filled circle represents the time-average 
generation rate of boron interstitials for the interval from s to 600 s 



It is worth noting that for simulation of boron redistribution during an- 
nealing of ion-implanted layers the mechanism of transient enhanced diffusion 
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due to formation, migration and dissociation of the "boron atom - silicon self- 
interstitial" pairs is usually used. However, we can neglect this mechanism 
of boron diffusion in simulating the experimental data of [6] due to the very 
low temperature of annealing and, hence, due to the very small thermal bud- 
get. Indeed, the boron diffusivity D for a temperature of 600 °C calculated 
according to the results of [23] is equal to 1.78xl0 -13 /im 2 s _1 . Then the 
average length of boron redistributions due to the pair diffusion mechanism 
L = ^/DTprC^/Cf 1 equals 1.4 x 10 3 /im (1.4 nm) for annealing duration of 
1200 s and for concentration of nonequilibrium silicon self-interstitials C AI is 
10000 times higher than the thermally equilibrium value of this species Cf 1 . 
Therefore, it is possible to neglect the traditional transient enhanced diffusion 
Really, one can see from Figs. [TJ and [2] that there is no boron redistribution 
in the high-concentration region of impurity. 

It is interesting to compare the derived value Iai = 24 nm for B implanted 
in crystalline silicon with the values of the average migration length of boron 
interstitials obtained for the case of B implantation in the layer of silicon 
preamorphized by Ge ions. For example, in [8] <100> oriented, Czochralski 
grown silicon substrates of the n-type of conductivity and resistivity of about 
10-25 Qm, were used for such experiments. These substrates were implanted 
with Ge ions at 32 keV to a dose of lxlO 15 cm -2 that provides amorphization 
to a depth of ~ 55 nm, as determined by Rutherford backscattering spectrom- 
etry (RBS) and cross-sectional transmission electron microscopy (XTEM). 
Boron was subsequently implanted at 500 eV to a dose of 2xl0 15 cm -2 at 
°tilt and °twist. A 1.5 MeV He beam at 45° to the sample were used for RBS. 
After implantation, samples received isochronal rapid thermal annealing in 
the N2 ambient at a temperature of 850 °C for 60 s. Atomic profiles of B 
were measured by SIMS. In [TTJ modeling of the doping process described 
above was carried out, and the average length of boron interstitials Iai = 12 
nm was obtained. In [9] boron implantation with an energy of 1 keV and 
a dose of 1.5 xlO 15 cm~ 2 was carried out in (100) oriented silicon substrates 
with the n-type of conductivity. For suppression of the transient enhanced 
diffusion a near surface layer of silicon was preamorphized by implantation 
of Ge ions with an energy of 15 eV. After boron implantation, the thermal 
annealing in the nitrogen ambient was carried out for 60 s at a temperature of 
800 °C. Calculation of the redistribution of ion-implanted boron for this ex- 
perimental environment performed in [13J has shown a good agreement with 
the experimental data for the average migration length of boron interstitials 
Iai = 11 nm. 
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It is clear from the presented results that the average migration length of 
interstitial boron atoms in crystalline silicon is significantly larger than in the 
silicon preamorphized by Ge ions and recrystallized due to the solid phase 
epitaxy at the initial stage of annealing. The obtained result confirms the 
strong interaction of migrating boron interstitials with Ge atoms previously 
implanted in the silicon layer. This interaction provides approximately a 
quadruple reduction of the average lifetime of boron interstitials and the 
relevant decrease in their migration length. 

6 Conclusions 

Simulation of impurity redistribution for boron implanted in crystalline sil- 
icon during annealing at a temperature of 600 °C for 600 and 1200 s has 
been carried out. The results obtained extend and deepen the investigation 
described in [15] where boron diffusion for 10-s annealing has been simu- 
lated. The boron concentration profiles calculated within the framework of 
the model of the long-range migration of impurity interstitials agree well with 
the experimental data for low-temperature thermal treatment of crystalline 
silicon. It is interesting to note that the same value of the average migra- 
tion length of nonequilibrium boron interstitials Iai = 24 nm can be used 
to explain, within the experimental error, the formation of extended "tails" 
in the low-concentration region of boron concentration profile for annealing 
duration of 10, 600, and 1200 s as well as for explaining the fast interstitial 
diffusion during ion implantation It was also shown that approximately 1.5 
% and 2.29 % of the implanted atoms occupied temporally the interstitial 
position during annealing for 600 and 1200 s, respectively, participated in the 
fast diffusion, and then became immobile again, took up the substitutional 
position or formed complexes with defects of the crystalline lattice. It follows 
from the obtained values of the relaxation time for generation of nonequi- 
librium impurity interstitials that the main fraction of boron interstitials is 
generated during annealing duration shorter than 30 s. 

A comparison of the present calculations with the earlier published data 
shows that the obtained value of Iai is approximately two times greater 
than the average migration length of boron interstitials in the silicon lay- 
ers preamorphized due to Ge ion implantation and recrystallized by means 
of solid phase epitaxy at the initial stage of annealing. It means that the 
migrating boron interstitials interact intensively with Ge atoms, which have 



13 



been used for the formation of a near-surface amorphous layer before boron 
implantation. 
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